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ABSTRACT 

This  paper  presents  a  reliable  control  design  technique 
for  linear,  time-invariant,  multi-input  multi-output 
(MIMO)  systems  with  degraded  actuators  and  sensors. 
The  degradation  defined  in  this  paper  ranges  from 
normal  operational  conditions  to  complete  failure  of 
actuators  and  sensors.  We  derive  linear  matrix  inequa¬ 
lity  (LMI)  conditions  ensuring  robust  stability  of  the 
system  using  static  state  feedback.  The  potential  of  the 
proposed  technique  has  been  demonstrated  by  an  ex¬ 
ample  of  three  coupled  inverted  pendulums. 

Keywords — linear  matrix  inequality;  static  state 
feedback;  MIMO;  degradation. 

I.  Introduction 

Conventional  feedback  control  designs  for  a  MIMO 
system  may  result  in  instability  in  the  event  of  the 
degradation  of  actuators  and  sensors,  even  though  it 
may  be  possible  to  control  the  system  using  only  the 
surviving  functions  of  actuator  and  sensor.  It  is  worth 
noting  that  this  condition  can  occur  even  if  the  open 
loop  system  is  stable.  It  is  therefore  of  interest  to 
develop  feedback  control  designs  which  guarantee 
robust  stability  despite  degradation  of  actuators  and 
sensors. 

There  are  relatively  few  methods  for  design  of  reliable 
control.  Veillette  et  al  [1]  develop  observer-based 
reliable  centralized  and  decentralized  control  systems 
provided  the  failure  of  actuators  and  sensors  occur 


only  within  a  prespecified  subset  of  the  control 
components.  In  addition,  the  actuators  and  sensors  are 
subject  to  either  failure  or  normal  operation  which 
neglects  the  possibility  of  partial  degradation.  Fujita 
and  Shimemura  [2]  provide  necessary  and  sufficient 
conditions  by  introducing  a  class  of  {/-matrices  for 
examining  system  stability  against  arbitrary  feedback- 
loop  failures  in  multivariable  control  systems  with  a 
stable  plant  and  controller.  This  condition  requires  an 
exponential  growth  in  computation  as  a  function  of 
plant  dimension.  Partial  failures  are  not  considered. 
For  a  given  stabilizing  controller  of  a  plant, 
Vidyasagar  and  Viswanadham  [3]  studied  the  problem 
of  designing  the  second  controller  such  that  either 
controller  acting  alone  stabilizes  the  plant  and  both 
acting  together  also  stabilize  the  plant.  A  computing 
method  for  the  second  stabilizing  controller  which 
involves  stable  coprime  factorizations  using  the  plant 
and  the  first  controller  data  are  established  and  may 
result  in  the  second  controller  of  high  order.  Cho  et  al 
[4]  used  the  same  methodology  as  [3]  except  that  the 
second  stabilizing  controller  is  of  adaptive  controller 
which  again  ends  up  with  high  order. 

In  this  paper,  we  design  a  static  state  feedback  control 
system  to  tolerate  the  degradation  ranging  from 
normal  operational  condition  to  complete  failure  of 
actuator  and  sensor.  The  degradation  is  modeled  as  a 
multiplicative  uncertainty  at  the  plant  input  or  output. 
We  accomplish  this  via  diagonal  weighting  and  norm 
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bounded  matrices.  To  simplify  the  derivation  we  will 
first  develop  and  prove  LMI  conditions  which  tolerate 
sensor  or  actuator  degradation  independently.  Our 
approach  is  extended  to  the  simultaneous  degradation 
actuators  and  sensors  in  the  system. 

II.  Problem  Formulation 


(3)  0  <  (p,!'  <  1  indicates  the  rth  actuator  may  fail 
partially. 


The  purpose  of  this  section  is  to  define  the  framework 
on  which  our  approach  of  reliable  control  systems  is 
based.  The  degradation  for  actuators  and  sensors  will 
be  considered  separately  to  simplify  the  derivation. 


A.  Actuator  case: 

Actuator  case  represents  the  degradation  of 
actuators  and  is  shown  as  follows, 


x  =  Ax  +  Bu 
u  =  uk  +  pu 


<lu 

r 

Pu 


(PuUk 


(2.1) 


where  xeR”  denotes  the  state  of  the  system,  uk  e  Rm  is 
the  output  of  the  control  gain  shown  in  Figure  1.  u  e 
R”‘  is  the  input  to  the  plant.  ip«  is  the  uncertainty 
weighting  of  the  control  input  to  the  plant.  Note  that 
pu  and  qu  represent  the  output  and  input  vectors  of 
plant  input  perturbations.  pur  and  qJ  correspond  to  the 
rth  element  of  vector  pu  and  qu  respectively. 


Figure  1  Closed  Loop  System  for  Actuator  Case 

B.  Sensor  case: 

Sensor  case  represents  the  degradation  of  sensors 
and  is  written  as  follows, 

x  =  Ax  +  Bu 

,  (2.2) 
xs  =  x  +  ps  v  ’ 

qs  =  Vs* 

pls  =  Slsqls,  |^(0|<1. 

where  y/s  is  the  uncertainty  weighting  of  the  sensor  and 
xs  denotes  sensor  output.  The  definition  for  ps,  qs,  y/s, 
and  8s  are  defined  in  (2.3)  in  a  similar  way  as  in  the 
actuator  case  and  can  be  seen  clearly  from  Figure  2 
where  it  shows  the  closed-loop  system  for  sensor 
degradation. 

p$  =  [psX  ....  psl ...  psn  ]T,  pJ^R 

qs  =  [  qs1  ....  qJ  ...  qs"  ]T,  qJe  R  (2.3) 


Pu  =  [puX  ....  Pu  . 

n  m  1 7 

■■  Pu  J  > 

pu  e  R 

K 

II 

K 

R 

■'S 

n  m  R 

.  qu  J  , 

qur  e  R 

q)u  =  diag  (  q>u  ... 

.  (pj ...  (pum). 

(pj  e  R 

5ur(t)  is  the  real  valued  unknown  bounded  uncertainty 
and  it  is  assumed  to  be  time-varying  scalar  quantity. 
The  control  signal  u(t)  can  be  readily  computed  as 


8s  =  diag  (A1 ....  8sl ...  8"  ),  Ss  e  R 

y/s  =  diag  (  y/sx ....  y/sl ...  y/A),  y/J  e  R 

The  degraded  sensor  signal  xs{t)  can  be  readily 
represented  as 

xs(t)  =  ( /  +  8s(t)  y/s)x(t) 


U(t)  =  (  I  +  8,(t)(p,i)l/k(t) 

where  the  / is  identity  matrix  and  <X  =  diag  ( 8u 1 . . . .  8/' ... 
8r  ),  81  e  R.  Since  <Xf(t)|  <  1,  i.e.  -1  <  8/(t)  <  1, 

the  degradation  is  modeled  by  a  priori  weighting  cpu  as 
follows. 

(1)  cpu  =  1  indicates  the  Ah  actuator  may  fail  fully. 

(2)  (pj  =  0  indicates  there  is  no  possibility  that  the  Ah 
actuator  fails  partially  or  fully. 


where  the  possibility  of  full  or  partial  failure  of  sensors 
is  defined  in  the  similar  way  as  in  the  actuator  case 
using  variables  Ss  and  y/s  in  place  of  8,  and  (pu, 
respectively. 


Figure  2  Closed  Loop  System  for  Sensor  Case 


@  IJTSRD  |  AvailableOnline@www.ijtsrd.comj  Volume- 1|  Issue -5 


Page:91 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 


Actuator  and  sensor  case 

Actuator  and  sensor  case  is  the  case  where  actuator 
and  sensor  are  simultaneously  degraded  shown  as 
follows, 

x  =  Ax  +  Bu 
u  =  uk  +  pu 
xs  =  X  +  ps 

qu  =  (puuk  (2.3) 

qs  =  V  sx 

p:  =  k<,  |*;co|si 

P\  =  8lq[,  1^(01  <  i. 


Figure  3  Closed-loop  System  for 
Actuator  and  Sensor  Case 

The  matrices  A  and  B  shown  in  (2.1),  (2.2),  and  (2.3) 
are  assumed  to  be  real  constant  matrices  with 
appropriate  dimensions.  (p»  and  i/a  are  given  a  priori 
and  represent  the  uncertainty  weighting  used  to  model 
the  degradation  for  actuators  and  sensors.  We  define 

A u  :=  diag  (5C (f))/=i 5 u(t)  e  R 

As  :=  diag  (8sl(t))i=i,...,n,  8 s‘(t)  e  R 

where  ||A«||X<  1  and  ||Aj||®<  1. 

The  control  problem  is  to  find,  if  it  exists,  static  state 
feedback  that  stabilizes  the  system  described  in  (2.1), 
(2.2),  or  (2.3). 

III.  Controller  Synthesis 

In  this  section  we  present  the  main  results.  Lyapunov 
stability  theory  is  used  to  establish  stability  criteria. 
Although  there  is  no  trivial  method  to  establish  a 
Lyapunov  function,  by  experience,  quadratic 
Lyapunov  functions,  i.e.  V(E)  =  c'Lc,  have  been 
proven  to  be  efficient  and  easily  implemented  [7].  We 
consider  static  state  feedback,  i.e.  uiAjt)  =  Kx(t)  for  the 
actuator  case,  u(t)  =  Kxs(t)  for  sensor  case  and  uk(t)  = 
Kxs(t )  for  combination  of  actuator  and  sensor  case. 


Theorem  1.  Consider  the  closed-loop  system 
x  —  (A  +  BK)x  +  Bpu 

q  u  =  <PuK x  t3-1) 

pru  =  5Wu>  k:(o|<i. 


If  there  exist  matrices  Q,  Y,  and  M  satisfying 
Q=QT>  0 

~QAT  +  AQ+BY+YtBt  +BMBt  Y'ipl 
<PJ  ~M 


<0 


(3-2) 


for  a  given  weighting  matrix  cpu,  then  the  following 
statements  are  equivalent, 

(1)  The  closed-loop  system  (3.1)  is  asymptotically 
stable. 

(2)  (A  +  BK)  is  Hurwitz  for  K  =  YQ~l. 

Proof.  See  Appendix  A  for  proof  and  notation. 

Remark  1.  The  inequality  (3.2)  is  convex  in  the  matrix 
variables  Q,  Y,  and  M.  Thus,  (3.2)  can  be  readily 
solved  by  convex  optimization  techniques.  We  write  a 
feasibility  problem 

Find  Q,  Y,  and  M 

Subject  to  Q  =  QT  >  0,  M>  0,  and  (3.2). 

Theorem  2.  Consider  the  closed-loop  system 
x  =  (A  +  BK)x  +  BKps 


qs  =  v  sx 

i  ?  i  i  L  i 


(3-3) 


Ps  =  8  sqs,  |8  s  (0|  -  1- 

If  there  exist  matrices  P,  W,  S  and  A  satisfying 
P  =  PT>  0,  S  =  A7'>0,  A  >  0, 

-PAP  +  WrS~lW<0  (3.4) 


PAT  +  AP+WTBT  +  BW  +  BSBt  P\\i  ts 


V/ 


- 


<  0 


(3.5) 


for  a  given  weighting  matrix  y/s,  then  the  following 
statements  are  equivalent. 

(1)  The  closed-loop  system  (3.3)  is  asymptotically 
stable. 

(2)  (A  +  BK)  is  Hurwitz  for  K  =  WP~l. 

Proof.  See  proof  for  Theorem  3. 
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Remark  2.  The  matrix  inequalities  (3.4)  and  (3.5)  are 
inversely  coupled  through  the  matrices  S  and  A. 
Moreover  the  inequality  (3.4)  is  not  an  LMI.  Thus,  to 
solve  theorem  2  simultaneously  through  the  use  of 
convex  optimization  is  difficult.  However,  we  may 
still  solve  the  problem  using  of  sequential  method  that 
we  solve  for  (3.5)  first  and  then  verily  (3.4). 

Theorem  3.  Consider  the  following  closed-loop 
system 

x  =  (A  +  BK)x  +  BKps  +  Bpu 
9u  =  (PuK(x  +  ps) 

(3.6) 

9,  =  ¥  sx 

p:  =  s’x„  K'(o|si 
p',  =  sw„  K«)|si. 

If  there  exist  matrices  P,  W,  Su,  Sx,  Ss,  A,  and  T 
satisfying 

P=PT>0,  su  =  s:> 0,  SX=S*> 0,  s,  =  sf> 0, 

A>o,r>o 


IV.  Numerical  Example 

Consider  the  system  in  Figure  3  consisting  of  three 
coupled  inverted  pendulums  of  point  masses  rrn,  and 
length  li.  The  pendulums  interact  via  three  springs  and 
three  dampers  of  stiffness  ky  and  damping  by;  i,  j  = 
1,2,3,  and  zfy.  The  distances  from  attached  point  of 
springs  and  dampers  to  the  platfonn  baseline  are  au 
The  system  data  are  shown  in  Table  1.  Using  this 
system,  we  will  demonstrate  several  examples  for 
actuator  and  sensor  degradation.  The  system  dynamics 
for  pendulums  are  written  in  the  general  form, 


-A  +  BtS~1B  <  0 

(3.6) 

-  pap  +  ivT(Sx'  +  s;'  +  ®rw<  o 

(3.7) 

f  pat  +  AP+WTBT  +  BW] 

-  \  wT 

K+bsxbt+su  ) 

W  -(®rS,0)r+®r)' 

p\ u 

(3.8) 

0 

<0 

\| )  ,P  0 

-  A-1 

where  chr=  (Pu7r(pu ,  for  a  given  weighting  matrix  (pu 
and  \|/s,  then  the  following  statements  are  equivalent, 

(1)  The  closed-loop  system  (3.6)  is  asymptotically 
stable. 

(2)  (A  +  BK)  is  Hurwitz  for  K  =  WP'1. 

Proof.  See  Appendix  B  for  proof. 

Remark  3.  Normally  speaking,  the  matrix  inequalities 
(3.6)-(3.8)  should  be  solved  simultaneously.  However, 
it  is  not  possible  to  solve  using  convex  programming 
due  to  inversely  coupled  terms,  Ss,  Su,  Sx,  and  A. 
Thus,  we  suggest  for  some  given  Or  we  solve  (3.8)  for 
P,  W,  Sx  ,  Ss,  and  Su  and  then  verify  (3.6)  and  (3.7). 


Figure  3  Three  Inverted  Pendulums  Systems 
Table  1 


rrn 

m2 

m3 

bn 

bi3 

b23 

ki2 

ki3 

k23 

1 

1 

1 

1 

.51 

1 

1 

.51 

1 

ai-a2=a3= 0.5,  g=10,  li=  1,  /^=  1 .2,  /.?=  1.1. 


Initial  conditions: 

ex  =-0.1,  d\  =0.11,  02  =0.34,  02  =-0.23,  <?3  =-0.29,  <?3  =-1.33. 

We  will  first  demonstrate  actuator  degradation.  The 
weighting  is  chosen  ^;,=diag(l,  0.1,  0.28)  which 
represents  the  possibility  of  degradation  for  each 
actuator.  For  instance,  actuator  for  rrn  may  fail  during 
its  operation,  actuators  for  mi  and  mi  are  subject  to 
10%  and  28%  variation  of  its  nominal  operation  signal. 
We  will  demonstrate  nominal  operation  of  actuators, 
i.e.  (fhi  =  0,  at  time  t  <  0.5.  Then,  we  will  show  the  case 
where  the  actuator  for  m\  is  subject  to  fail  and  the 
actuator  for  mi  and  mi  are  partial  failure.  Detail 
conditions  are  shown  in  Table  2. 


Table  2 


time 

m  i 

m2 

m  3 

0  <t  <  0.5 

Si 

S2 

S3 
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0.5  <t  <  2.5 

failure 

0.9*S2 

0.72*Ss 

2.5  <  f  <  6 

2*Si 

1.1*S2 

1.28*S3 

chosen  y/s  =  diag(0.4,  0.5,  0.4,  1,  0.5,  1)  which 
represents  the  possibility  of 


where  S«  =  the  nominal  operation  control  signal  of 
actuator  associated  with  each  mass.  The  computed 
gain  followed  by  Theorem  1  is 


K  = 


-44.6423  -11.3417  2.7025  -1.3565  1.2500  -0.7005 

-362.5824  -93.3584  -2.9357  -17.3282  -4.3807  -10.9399 

-25.4243  -7.1539  4.5474  -1.6523  -35.8191  -10.0032 


eigenvalues  of  (A+BK)  =  [-45.2643  -3.6514 

1.13 16+1 .6097i  -2.13 06+2 . 1 3  94i] .  The  simulation 
results  are  shown  in  Figure  4  and  5  for  plant  outputs 
and  control  signals.  It  is  noticed  that  the  control  signal 
for  mi  decays  relatively  fast  at  the  beginning  of 
operation  and  then  is  close  to  zero  as  shown  in  Figure 
4.  Thus,  even  though  the  actuator  #1  fails  at  time  t  > 
0.5  sec,  it  has  very  little  effect  to  the  system.  Since  the 
actuator  signals  for  m  i  varies  from  twice  of  its  nominal 
value  to  zero  (which  is  used  to  represent  the  failure  of 
the  actuator),  one  possible  way  to  accommodate  such 
a  large  variation  is  simply  to  make  the  resultant  signal 
as  small  as  possible  such  that  it  has  the  least  effect  on 


Figure  5  Control  Signals 

degradation  for  each  sensor.  The  velocity  sensor  of  mi 
is  subject  ++to  50%  variation,  while  the  position 
sensor  has  40%  variation.  The  velocity  sensors  of  mi 
and  tm  are  subject  to  the  possibility  of  failure,  while 
the  corresponding  position  sensors  have  40%  and  50% 
variations.  We  will  demonstrate  the  extreme  case, 
which  is  shown  in  Table  3.  The  computed  gain 
followed  by  Theorem  2  is 


the  system.  The  jumps  shown  in  Figure  5  are  due  to 
the  switching  at  time  t  =  0.5  sec  instantaneously.  If 
there  are  two  actuators  in  the  system  subject  to  the 
possibilities  of  failure,  e.g.  ^»=diag(l,  1,  0),  we  are  not 
able  to  find  feasible  solutions  for  the  Theorem  1.  In 


K  = 


-29.7241 

0.5669 

1.2812 


-8.2556  1.0801 

-0.9255  -29.5401 
-0.2544  1.2151 


-1.0596  1.2811  0.3397 

-8.0746  0.9343  0.9577 

-0.9838  -29.3428  8.3852 


eigenvalues  of  ( A+BK)  =  [-6.8259  -2.4785  -5.8912  - 
3.0449  -3.0952  -6.2916].  The  simulation  results  are 


other  words,  it  is  not  possible  to  maintain  stability  of 
the  system  when  two  actuators  failures  in  the  system. 


Figure  4  Outputs  of  Three  Pendulums 
For  the  sensor  degradation  case,  the  weighting  is 


shown  in  Figure  6  and  7.  The  sudden  jumps  in  Figure 
7  show  the  instantaneous  switching  at  time  t  =  0.5  sec. 
The  following  results  are  observed  from  simulation. 
The  system  will  oscillate  if  we  fail  three  velocity 
sensors.  Since  the  damping  control  signals  given  by 
velocity  sensor  are  zero,  the  occurrence  of  oscillation 
is  not  unexpected. 


Table  3 


time 

mi 

mi 

m3 

Pos 

Vel 

Pos 

Vel 

Pos 

Vel 

0<  t 

1.4* 

1.5* 

1.4* 

2* 

1.5* 

2* 

<  .5 

Sip 

Slv 

S2p 

S2v 

S3p 

S3v 

0.5< 

0.6* 

0.5* 

0.6* 

F 

0.5* 

F 

t<  5 

Sip 

Slv 

S2p 

S3p 

Pos  =  position  sensor,  Vel  =  velocity  sensor,  S  with  all 
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subscripts  =  the  nominal  operational  signal  associated 
with  each  sensor  (subscript  p  =  position  sensor, 
subscript  v  =  velocity  sensor).  F  =  failure. 


Figure  6  Outputs  of  Three  Pendulums 


Proof. 

We  consider  the  quadratic  Lyapunov  function  V(x)= 
xTLx.  The  consideration  of  perturbation  will  be 
incorporated  in  the  derivative  of  Lyapunov  function 
for  stabilization  of  overall  system.  We  have 


V(x)  +  YJ  Lr(q'„Tq'u  -  pi! p'u)  <  0  (A.  1) 

r—\ 

for  some  X,->0.  Substituting  (2.1)  into  (A.l)  and 
rearranging  the  expression,  we  obtain 


x 

Pu 


'  AtL+ Li + LBK + ( LBK ) r 

LB 

X 

I+kYM.k 

V 

B‘L 

-A 

_LU_ 

(A.2) 


where  the  A=diag(ki,  Xi, ...,  Xm).  (A.2)  is  equivalent  to 


Figure  7  Control  Signals 


V.  Conclusions 

We  ensure  the  stability  of  a  MIMO  system  by 
establishing  LMI  conditions  when  the  system  is 
subject  to  degradation  of  actuators  and  sensors.  The 
theory  developed  in  Section  3  is  demonstrated  by  three 
inverted  coupled  pendulums,  which  show  the  system 
can  be  stabilized. 
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Appendix  A 


AT  L  +  LA  +  LBK  +  (LBK) 
V+KT<pTuX<puK 


t\ 


btl 


LB 

-A 


<  0 


(A. 3) 


Expanding  by  Schur  complement  and  letting  Q  =  L~l, 
Y  =  KQ,  and  M  =  A'1,  we  may  have 

~QAT  +AQ+BY+YtBt  +BMBt 

<pj  -M  J<0 

This  completes  the  proof. 

Appendix  B 

Proof. 

We  need  the  following  lemma  used  in  the  proof. 

Lemma  B  ([8])  For  any  matrices  M  and  N  with 
appropriate  dimensions,  we  have 

MtN  +  NtM<  MTSN+NTS~lM 

for  any  positive-definite  symmetric  matrices  S. 

We  consider  the  quadratic  Lyapunov  function  V(x)= 
xTLx.  The  system  perturbation  will  be  incorporated  in 
the  derivative  of  Lyapunov  function  for  stabilization 
of  overall  system.  We  have 


r=l 


m 

-prsTp's)+ Yj?M'«qr„ -PuPru)  <  0 
r= 1 


(B.l) 
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for  some  U>0  and  6>0.  Substituting  (3.6)  into  (B.l) 
and  rearranging  the  expression,  we  obtain 


xT (A[ L  +  LAt  +  KylTy^K  +  yirA\|/l)x  N 
+  xT  (LB  K  +  KT<pTuryiiK)ps  +  x  T  LBp  u 
+  p](KTBT  L  +  A>,(r(pi(A)x  +  pIBtLx 
,+  pi (-A  +  KT  yif  q  uK)  p  s  -  plTpu  , 


(B.2) 


where  the  A=diag(ki,  ki, ...,  km)  and  T=diag(t i ,  £2, 
£,m).  Use  of  Lemma  B  we  have 

xT LBp u  +  pi Bt Lx  <  xTLSIILx  +  pTuBTS;'Bpu 
xTLBKps  +  pi KtBt Lx  <  xTLBSxBTLx  +  pTs  KrS~x'KPs 
xTKT®TrKps  +  pTs  KT®TrKx  <  xTKT®TrSs®rKx  +  pTs  KTS~lKps 

where  or  =  cp,£rcp„ .  Thus, 


(B.2)  < 


'  J  ATkL  +  LAk  +KT(®r  +®TrSs®r)K  +  ylAy^ 
x{  +l(Su  +  bsxbt)l  Jx 

,  +  P:(-A  +  Kr(Sd  +S-'+®Y)K)Ps+PTu{-T  +  BTS;'B)Puy 


<0 


Therefore,  we  require 

^IL  +  LAk  +  /:r((l)r  +  r)^l  <  0  (B.4) 

l+V.rA  vs+L(Su+BSxBt)L  )< 


-A  +  KT(S;l+S;l+d>r)K<0  (B.5) 

-r  +  Rr5,;1R<0  (B.6) 


Multiplying  P  =  Z"1  on  both  side  of  (B.4)  and  applying 
Schur  complement,  we  obtain 
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-A'1 

(B-7) 


where  W=KP.  Multiplying  P  on  both  side  of  (B.5)  we 
have 


-pap+wt(s;'  +s;'  +®rw<o  (B.8) 

This  completes  the  proof. 
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